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Straintronic logic is a promising platform for beyond Moore’s law computing. Using Bennett clocking mech¬ 
anism, information can propagate through an array of strain-mediated multiferroic nanomagnets exploiting 
the dipolar coupling between the magnets without having to physically interconnect them. Here we perform 
a critical analysis of switching failures, i.e., error in information propagation due to thermal fluctuations 
through a chain of such straintronic devices. We solved stochastic Landau-Lifshitz-Gilbert equation consider¬ 
ing room-temperature thermal perturbations and show that magnetization switching may fail due to inherent 
magnetization dynamics accompanied by thermally broadened switching delay distribution. Avenues available 
to circumvent such issue are proposed. 


Multiferroic devices consisting of a piezoelectric layer 
strain-coupled to a magnetostrictive nanomagnet, hold 
profound promise to replace traditional transistors for 
our future information processing paradigm. These de¬ 
vices work according to the principle of converse magne¬ 
toelectric effecti.e., when a voltage is applied across 
the device, the piezoelectric layer gets strained and the 
strain is elastically transferred to the magnetostrictive 
layer rotating its magnetization (see Fig. [1^). With 
appropriate choice of materials, such devices dissipate 
a minuscule amount of energy of attojoule in sub¬ 
nanosecond switching delay at room-temperature.— This 
study has opened up a field called straintronics^ and 
experimental efforts to demonstrate such electric-field in¬ 
duced magnetization switching are considerably emerg¬ 
ing jSi^— 

Information processing using Bennett clocking mecha- 
nisn>ii is an attractive platform for building logic using 
nanomagnets and dipolar coupling between them. This 
facilitates avoiding physical interconnects and thus elim¬ 
inating the energy dissipation due to charging and dis¬ 
charging of interconnect capacitances. Figure [TJj depicts 
how a bit of information can be propagated unidirection- 
ally through a chain of energy-efficient stress-mediated 
multiferroic devices^ i^^’^^ rather than using highly energy 
consuming magnetic heldji^^i^ Contrary to the steady- 
state analysisfi^ the investigation of magnetization dy- 
namicsii^i^ has proved to be crucial for achieving sub¬ 
nanosecond switching speed making the straintronic logic 
competitive with traditional charge-based computing.— 
Recent experiments have demonstrated defects and er¬ 
rors in nanomagnetic logic circuits jdi 

Here, we study the source of switching failures in 
straintronic logic due to thermal fluctuations during the 
propagation of a bit of information. We solved stochastic 
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Landau-Lifshitz-Gilbert (LLG) equation^d2r.2£ of magne¬ 
tization dynamics to understand the critical issues be¬ 
hind the switching failures. As such we would assume 
that the magnetization of the 2nd nanomagnet in Fig.jTjD 
would always switch successfully to the desired state as 
shown in the last row due to the dipole coupling from the 
1st nanomagnet. However, the analysis presented here 
demonstrates that magnetization’s slight excursion out 
of magnet’s plane accompanied by the thermal fluctua¬ 
tions can eventually make magnetization backtracking to 
the wrong direction. This would produce error in prop¬ 
agating a bit of information. Making an approximation 
by not taking into account the out-of-plane excursion of 
magnetization2ii2^ would not be able to comprehend such 
critical reasoning behind switching failures. Noting that 
it requires a very small bit error rate (< 10“^) for com¬ 
puting purposes, we further suggest an way to tackle such 
issue. 

We model the magnetostrictive nanomagnet in the 
shape of an elliptical cylinder; its cross-section lies on the 
y-z plane, the major axis points along the z-direction, 
and the minor axis along the y-direction (see Fig. [T^). 
Any deflection of magnetization out of magnet’s plane {y- 
z plane, (j) = ±90°) is termed as out-of-plane excursion. 
The dimensions of the major axis, the minor axis, and 
the thickness are a, b, and I, respectively {a > b > 1). So 
the volume is H = (7r/4)a6Z. We will consider the switch¬ 
ing of nanomagnet-2 and the subscript of any parameter 
will point to the corresponding nanomagnet (1 to 4, see 
Fig. [It). 

We solve the stochastic Landau-Lifshitz-Gilbert (LLG) 
equatioiii^^— in the presence of thermal fluctuations (de¬ 
tails are provided in the supplementary material^^) and 
derive the following coupled equations for the dynamics 
of 02 and (j) 2 - 
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FIG. 1. Bennett clocking mechanism for unidirectional information propagation in straintronic logic, (a) A voltage-controlled 
strain-mediated multiferroic device and axis assignment. Magnetization is bistable along the ± 2 ;-axis, which stores a bit 
of information 0 or 1. (b) Unidirectional information propagation through a horizontal chain of straintronic devices. The 
nanomagnets are stressed separately using different voltage sources. Note that the dipolar coupling between the neighboring 
nanomagnets is bidirectional and hence we need to impose the unidirectionality in time (using a 3-phase clocking scheme to 
apply stress on the nanomagnets subsequently) to propagate a bit of information through the chain. The magnetization of the 
1st nanomagnet is flipped, and the 2nd and 3rd nanomagnets are stressed to align their magnetizations along the hard-axis. 
Then stress is released/reversed on the 2nd nanomagnet to relax its magnetization along its desired state. 
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a is the phenomenological damping parameter, 7 is the 
gyromagnetic ratio for electrons. My = fioMsil, Ms 
is the saturation magnetization, Nd-mm is the compo¬ 
nent of demagnetization factor along m-direction, which 
depends on the nanomagnet’s dimensions ^^*2^ (3/2)As 
is the magnetostrictive coefficient of the single-domain 


magnetostrictive nanomagnet (T 2 is the stress on the 
nanomagnet -2 (note that the product of magnetostric¬ 
tive coefficient and stress needs to be negative in sign 
for stress-anisotropy to overcome the shape-anisotropy), 
At is the simulation time-step, G(o,i) is a Gaussian 
distribution with zero mean and unit variance^ k is 
the Boltzmann constant, T is temperature, Tdipoie,e 2 = 
{\ / silfl 92 '){dEdipole ,2 j ^^ 2 ) 7 Bdipole,(f >2 — 0 Prfzpo/e, 2 , 
and Edipoie ,2 is the dipole coupling energy from the neigh¬ 
boring nanomagnets 1 and 3. Note that in a very similar 
way the equations of dynamics for the other three nano¬ 
magnets can be derived. 

The magnetostrictive layer is considered to be made of 
polycrystalline Terfenol-D, which has the following ma¬ 
terial properties - Young’s modulus (Y): 80 GPa, satura¬ 
tion magnetization (Mg): 8x10® A/m, Gilbert’s damp¬ 
ing constant (a): 0 . 1 , and magnetostrictive coefficient 
((3/2)As): -l-90xl0~® (Refs. Isl. 12714^ . The dimensions 
of the nanomagnet are chosen as a = 100 nm, b — 90 nm, 
and I = 6 nm, ensuring the validity of single-domain as- 
sumption] ^®i®° The center-to-center distance between the 
nanomagnets is chosen as P = 120 nm. 

The piezoelectric layer is made of PMN-PT^ which 
has a dielectric constant of 1000 and the layer is assumed 
to be four times thicker than the magnetostrictive layer 
Assuming that maximum strain that can be generated in 
the piezoelectric layer is 500 ppm )®^0® it would require an 
electric field of ^^0.4 MV/m because ^31 = 13x10“^° m/V 
for PMN-PTj^ The stress generated in the Terfenol- 
D layer is the product of strain and Young’s modulus. 
Hence, 4.6 mVs of voltages would generate 20 MPa stress 
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Distribution of delay [9=90 ] under thermal fluctuations 



FIG. 2. Distribution of switching delay when magnetization of 
nanomagnet-2 reaches at 6 = 90° from 0 ~ 180° upon appli¬ 
cation of 20 MPa stress at 100 ps ramp period. A moderately 
large number (10,000) of simulations have been performed in 
the presence of room-temperature (300 K) thermal fluctua¬ 
tions to generate this distribution. This wide distribution is 
caused by the following two reasons: (1) thermal fluctuations 
make the initial orientation of magnetization a distribution, 
and (2) thermal kicks during the transition from 6 ~ 180° to 
0 = 90° make the time-period a distribution too. The mean 
and standard deviation of this distribution are 0.232 ns and 
0.056 ns, respectively. 

in the Terfenol-D layer. 

Figure [5] shows that upon application of stress, dif¬ 
ferent trajectories of magnetization of nanomagnet-2 
reaches 9 = 90° at variable times in the presence of ther¬ 
mal fluctuations. We take the distribution of initial orien¬ 
tation of magnetization due to thermal fluctuations into 
account. If we release/reverse the stress on nanomagnet- 
2 ahead of time, then magnetization may not be able to 
switch successfully rather it will backtrack to the same 
easy axis it started. This is exemplified in Fig. [3^. Mag¬ 
netization failed to switch to 6* ~ 0°, even it could not get 
past 9 ~ 90° since stress is kept constant for a short pe¬ 
riod of time of 100 ps. Out of 10,000 simulations, 16.52% 
switching failures were observed in this case. If we keep 
the stress constant longer for 200 ps, it would not nec¬ 
essarily result in reducing the failure rate of switching, 
which we will discuss later. Figure [5 }d depicts a case 
when switching failed; magnetization backtracked when 
stress was ramped down. We keep the stress constant for 
much longer time (600 ps) to observe whether sometimes 
switching still fails or not; Fig. |3t shows a case when 
switching fails. Note that magnetization keeps lingering 
around 9 = 90° since stress was kept constant for long 
time, but magnetization eventually backtracked towards 
9 = 180°. 

Figured^ depicts the non-trivial dependence of switch¬ 
ing failure rate with stress constant time. This can be 
explained by considering magnetization’s excursion out 
of magnet’s plane {y-z plane, (j) = ±90°) during switch¬ 
ing as shown in Fig. HJd. Magnetization deflects out-of- 
plane due to the torque exerted on it in the direc¬ 
tion and fast (non-adiabatic) ramp of stress. Although 
the dipole coupling from nanomagnet-1 facilitates switch¬ 
ing of magnetization of nanomagnet-2 towards 0 = 0°, 


the out-of-plane excursion of nanomagnet-2’s magneti¬ 
zation can hinder the switching. In Fig. |41 d, note that 
if magnetization resides in the quadrant (p £ (90°, 180°) 
[or (p £ (270°,360°)], the term Ss/iape.^a(‘/’z) as in the 
Eq. ([3]) becomes negative in sign, which would facilitate 
decreasing the value of 02 [see Eq. O] aiding magnetiza¬ 
tion rotation in the correct direction. But, if magnetiza¬ 
tion resides in the other two quadrants [cp £ (90°, 180°) 
or (p G (270°, 360°), termed as bad quadrants onwards], 
the term Bshape, 4 , 2 i 4 ’ 2 ) would be positive, which would 
force magnetization backtracking towards 0 = 180°. 
This inherent motion is generated particularly due to (p- 
dependence of potential energy, which is strong enough 
to affect the magnetization dynamics. This motion is also 
responsible for reducing the switching delay by a couple 
of orders in magnitude and bellied under the fact that 
^d—xx ^ ^d-yy-jd^d—zz {^d—xx ^d—yy Is higher than 
Nd-yy — Nd-zz by a couple of orders in magnitude). 

When magnetization reaches 0 = 90°, thermal fluc¬ 
tuations can scuttle magnetization in either side of the 
magnet’s plane {(p = ±90°). As explained earlier, if mag¬ 
netization happens to be in the had quadrants, magneti¬ 
zation would be prone to backtracking, causing a switch¬ 
ing failure. In Fig. |4^, when the stress constant time 
is low, any increase in stress constant time has two ef¬ 
fects that counter each other: (1) More trajectories fin¬ 
ish before stress is released (see Fig. W) so magnetiza¬ 
tion for those trajectories keep lingering around 0 = 90° 
and thermal fluctuations may cause mishap by scuttling 
magnetization in the had quadrants for (p. This will in¬ 
crease the switching failure rate. (2) More trajectories 
will be able to reach at 0 = 90° before we release the 
stress (see Fig. [2]), whereupon magnetization can possi¬ 
bly switch successfully towards 0 = 0°. This will decrease 
the switching failure rate. Due to these two counteract¬ 
ing effects, we do see a peak in the Fig. Now, if we 
keep increasing the stress constant time so that all the 
trajectories reach at 0 = 90° before stress is released, 
the switching failure rate saturates to '^7%. The reason 
is that thermal fluctuations scuttle magnetization in the 
bad quadrants for (p upon reaching 0 = 90° and dipole 
coupling from nanomagnet-1 cannot help much. This 
saturated switching failure rate of ~7% can be decreased 
a bit by optimizing parameters but it is not possible to 
decrease it to a very small value (< 10“"'^) for general- 
purpose computing purposes. Also, waiting until all tra¬ 
jectories reach at 0 = 90°, increases the switching time 
period. Hence, we propose the following approach. 

We realize that letting magnetization collapse on the 
magnet’s plane and thermal fluctuations scuttling mag¬ 
netization in the bad quadrants are causing the switch¬ 
ing failures. Provided we apply a sufficiently high stress 
with a sufficiently fast ramp rate, magnetization will not 
traverse into bad quadrants during switching while mag¬ 
netization passes through the x-y planei^ We can use a 
sensing circuitry to detect when magnetization reaches 
around 0 = 90° so that we can ramp down the stress 
thereafter. The sensing circuitry can be implemented by 
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FIG. 3. Magnetization of nanomagnet-2 backtracks to the same easy axis it started resulting in an error in propagating a bit 
of information, (a) Stress is ramped up from zero to 20 MPa in 100 ps, kept constant for 100 ps, and ramped down (and 
reversed) at the same rate as for ramp-up. Magnetization failed to even reach a.t 9 = 90° since stress was not active for a 
sufficient amount of time. Note the <()-dynamics that magnetization started close from an in-plane angle (j) = 270°(—90°), and 
traversed to another in-plane angle (f) = 90° due to 6-(j> coupled dynamics, (b) Stress is kept constant longer (for 200 ps) 
and ramp up/down times are same as for the part (a). Magnetization switching failure is still observed. However, this time 
magnetization was able to get past 9 = 90° towards 0 = 0°, but during the ramp-down phase, magnetization was subjected 
to a detrimental motion forcing magnetization to backtrack towards 9 = 180°. (c) Stress is kept constant even longer (for 600 
ps) and ramp up/down times are same as for the part (a). Magnetization switching failure is still observed. Magnetization is 
lingering around 9 = 90° since stress is kept constant for long time, however, as stress is brought down, magnetization started 
backtracking towards 9 = 180°. 
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Distribution of switching delay under thermal fluctuations 



FIG. 4. (a) Percentage switching failures versus stress constant time for 20 MPa stress and 100 ps ramp period. Failure 
rate increases with stress constant time initially but later decreases with further increasing stress constant time before getting 
saturated, (b) Out-of-plane {(j) ^ ±90°) excursion of magnetization during its dynamical motion. The magnetization deflects 
out-of-plane due to the torque exerted on it and fast (not adiabatic) ramp rate. When magnetization reaches around 9 — 90° 
{x-y plane), if it resides in the quadrant (90°,180°) [or (270°,360°)], it will aid magnetization’s motion towards 9 = 0°; 
otherwise, residing in the other two quadrants would be tantamount to switching failure, (c) Stress is released dynamically 
when magnetization reaches 9 — 90°. Magnetization always switches successfully in this methodology. The mean and standard 
deviation of this distribution are 0.459 ns and 0.067 ns, respectively. 


measuring the magnetoresistance in a magnetic tunnel 
junction (MTJ)i^i2ir— Basically, we need to get cali¬ 
brated on the magnetoresistance of the MTJ when mag¬ 
netization resides on the x-y plane. Comparing this 
known signal with the sensed signal of the MTJ, the 
stress can be ramped down. Figure|4]: shows the distribu¬ 
tion of switching delay considering such sensing circuitry. 
No switching failures were observed and the mean energy 
dissipation in the nanomagnets turns out to be ^^1.5 atto- 
joules at sub-nanosecond switching delay. Some tolerance 
is nonetheless required since the sensing circuitry cannot 
be perfect. We performed simulations to show that in¬ 


ternal dynamics works correctly as long as the stress is 
ramped down when magnetization’s orientation is in the 
interval 9 € (85°, 140°), i.e. it does not have to be exactly 
90°. This tolerance is due to the motion arising from the 
out-of-plane excursion of magnetization. 

In conclusion, we have performed a critical analysis of 
switching failures in energy-efficient straintronic logic us¬ 
ing Bennett clocking for computing purposes. It is shown 
that the switching failures are caused by the inherent 
magnetization dynamics particularly due to out-of-plane 
excursion of magnetization and thermal fluctuations dur¬ 
ing switching. We have proposed a novel approach to 
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circumvent such basic issue after a thorough analysis. 
Such methodology can be exploited for building logic 
gates and general-purpose computing purposes. Bennett 
clocking based architecture is regular in nature, so that 
different building blocks for computing purposes can be 
designed systematically. Such energy-efficient, fast, and 
non-volatile (that can lead to instant turn-on computer) 
computing methodology has profound promise of being 
the staple of our future information processing paradigm. 
Processors based on this paradigm may be suitable for 
applications that need to be run from energy harvested 
from the environment e.g., wireless sensor networks, med¬ 
ically implanted devices monitoring epileptic patient’s 
brain to warn an impending seizure. 
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Solution of stochastic Landau-Lifshitz-Gilbert (LEG) 
equation of magnetization dynamics in the presence of 
thermal fluctuations.- We adopt the standard spheri¬ 
cal coordinate system (see Fig. la in the main Letter) 
to solve the magnetization dynamics using stochastic 
Landau-Lifshitz-Gilbert (LLG) equation in the presence 
of room-temperature (300 K) thermal fluctuations^^ The 
potential energy of the nanomagnet -2 can be expressed 
as the sum of the anisotropies due to shape and stress, 
and due to dipolar coupling from the nanomagnets 1 and 

3, 


Gitotal,2 — Eshape,2 “f Estress,2 T Fldzpo/e,2 


= B2{(j)2) sin‘^92 -\- Edipole,2, 

(SI) 

where 


^ 2 ( 02 ) — Pshape,2i4^2) P Pstress,2i 

(S2a) 

Pshape,2{<t>2) = (Mo/2)Mg il[(A^d_yy — Nd- 

-zz) 

P {P^d—xx P^d—yy) 

005^02], 

(S2b) 

Pstress,2 — (3/2)As(T2il, 

(S2c) 


/^U o Q r 

Edipole,2 = Mg fl [cOs 02 COsdi COs 02 COsd 3 

sin9isin62{cos4>iCos4>2 — 2sin(j)isin(j)2) 
sin63sin62{cos4>3Cos<f>2 — 2sm(^3sm())2)], 

(S3) 

Ms is the saturation magnetization, Nd-mm is the 
component of demagnetization factor along m-direction, 
which depends on the nanomagnet’s dimensions,—*^ 
(3/2)As is the magnetostrictive coefficient of the single¬ 
domain magnetostrictive nanomagnet (J 2 is the stress 
on the nanomagnet-2, and R is the center-to-center dis¬ 
tance between the nanomagnets. Note that the product 
of magnetostrictive coefficient and stress needs to be neg¬ 
ative in sign for stress-anisotropy to overcome the shape- 
anisotropy. 


The magnetization M 2 of the nanomagnet-2 has a 
constant magnitude but a variable direction, so that 
we can represent it by a vector of unit norm rim ,2 = 
M 2 /IM 2 I = Br where Sr is the unit vector in the ra¬ 
dial direction in spherical coordinate system represented 
by {r,9,(j)). The effective held and torque acting on the 
magnetization due to gradient of potential landscape as 
in Eq. (ISTI) can be expressed as Hefr ,2 = —^Etotai ,2 = 
-{dEtotai,2/d92)eg - {l/sin92) {dEtotai, 2 /d(j) 2 )e^ and 
Te ,2 = rim, 2 X Heff, 2 , respectively. The thermal held 
and the corresponding torque acting on the magnetiza¬ 
tion can be written as Hth ,2 = eg -\- e,f, and 
Tth .2 = rim. 2 X Hth. 2 , respectively^*^ where 

Pg^ = Mv[hx ,2 cos 92 cos(j )2 + hy ^2 cos92sin(j)2 

— hz,2sin92], (S4a) 

P 02 = Mv[hy ^2 cos(j )2 — hx ,2 sm()) 2 ], (S4b) 

= = (S4c) 


a is the phenomenological damping parameter, 7 is the 
gyromagnetic ratio for electrons. My = iiqMsUI, At is 
the simulation time-step, G(o 7 ) is a Gaussian distribution 
with zero mean and unit variance^ k is the Boltzmann 
constant, and T is temperature. 

The magnetization dynamics under the two aforesaid 
torques is described by the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation^i^*21 as 




‘-111..^ 


dt 


- a rim ,2 X 


dt 


l7l . 


Tth,: 


(S5) 

After solving the LLG equation, we get the following cou¬ 
pled equations for the dynamics of 02 and 4)2 ■ 


(1 T Q! ) — [Bshape,<p2{'t’2)sin92 

at IVly 

2 g;.S 2 (02 '^dipole.,92 ^'^dipole,4>2 

+ {aPe^ + ^(^ 2 )]; (S6) 
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where 

/ , ^ dBshape,2{4>2) 

shape,<p 2 \Y 2 .) — ~ 

= (/io/2) MsVt{Nd-xx - Nd-yy)sin(2(f>2), 

(S8) 

'^dipole ,02 — /S^‘hl92)(,dEdipole,2/^^ 2 ^): B-Uci Tdipole ,^2 — 

dEdipoie, 2 /992- Note that in a very similar way the equa¬ 
tions of dynamics for the other three nanomagnets can 
be derived. 

The energy dissipation in the nanomagnet-2 due to 
Gilbert damping can be expressed as Ed ,2 = fj Pd, 2 {t) dt, 
where r is the switching delay and the instantaneous 
power dissipation can be calculated as 

While calculating internal energy dissipation, we sum 
up the energy dissipations in all the four nanomagnets 
but note that the dissipations in the nanomagnet-1 and 
nanomagnet-4 are quite negligible since they don’t rotate 
much and the dissipation in nanomagnet-3 is about half 
of that in nanomagnet-2 (since nanomagnet-2 switches 
a complete 180°, while the nanomagnet-3 switches only 
about 90°). We sum up this internal energy dissipa¬ 
tion with the energy dissipation in the external circuitry 


(which is miniscul o^^'^^ ) to get the total energy dissipa¬ 
tion. 

The performance metrics switching delay and en¬ 
ergy dissipation are determined by solving stochas¬ 
tic Landau-Lifshitz-Gilbert equation in the presence of 
room-temperature (300 K) thermal fluctuations. We as¬ 
sume a 20 MPa of stress and we reverse the stress during 
ramp-down phase to aid switching)^ We perform a mod¬ 
erately large number (10000) of simulations and we con¬ 
sider that magnetization’s initial orientation is not fixed 
rather it is fluctuating due to thermal agitations, so we 
have a distribution of magnetization’s initial orientation 
when magnetization starts switching. We determine the 
initial distributions of both polar angle 9 and azimuthal 
angle (j) at room-temperature by solving the stochastic 
LLG equation when no stress is active j^^i^^ When mag¬ 
netization of nanomagnet-2 reaches 9 < 5°, the switching 
is deemed to have completed. 
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